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A comprehensive analysis of the lepton flavor violating (LFV) decays Z — > is presented 

within the effective Lagrangian approach. Both the decoupling and nondecoupling scenarios are 
explored. The experimental constraints from li —> ijik^k and ii —» Ij^ as well as some relationships 
arising from the gauge invariance of the effective Lagrangian are used to put constraints on Z —> 
if l^. It is found that while current experimental data impose very strong constraints on Z ^ 
the channel Z — > still may be at the reach of the planned TESLA collider. 
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I. INTRODUCTION 

Recent neutrino experimental data, such as those coming from Super-Kamiokande Q, have shown evidences of 
atmospheric and solar neutrino oscillations. This class of effects point to physics beyond the standard model (SM) 
and have immediate consequences on some sectors of the theory. For instance, the conservation of lepton number 
and lepton flavor can not be taken for granted anymore, as in the SM with its massless neutrinos. Clearly, some 
lepton flavor violating (LFV) processes such as Z ^ P’ '*’) occur and be observable at the future 

particle colliders. The neutrino experimental data have thus renewed the interest in LFV transitions. Moreover, the 
prospect of the e+e“ TeV energy superconducting linear accelerator (TESLA) with its Giga-Z option opens up 


the possibility of studying at a near future some LFV Z boson decays, which might be at the reach of that collider 
[^. Currently, the best direct experimental bounds on the Z ^ if rates, obtained by the search at LEP-I, are Q 

BR(Z - 

e^p.^) < 1.7 X 10“®, 

(la) 

BR(Z - 

e^T^) < 9.8 X 10“®, 

(lb) 

BR(Z- 

4 pAr^) < 1.2 X 10“®, 

(Ic) 

whereas the expectations at TESLA are 



BR(Z - 

e^p^) < 2.0 X 10“®, 

(2a) 

BR (Z ^ 

T=^) < (1.3-6.25) X 10“®, 

(2b) 

BR (Z ^ pL^ 

T^) < (0.44-2.2) X 10“®. 

(2c) 


In order to disentangle the origin of any possible LFV effect, TESLA expectations must be confronted with the 
predictions of the diverse available models. Considerable work has been done along these lines, but here we will only 
refer to the most recent studies. For instance, the authors of Ref. reviewed diverse scenarios that enlarge the SM 
by just adding massive neutrinos. It turns out that, after considering the most recent experimental data for neutrino 
masses, one can have at most BR (Z If ^ 10“®"* for light neutrinos. On the other hand, Z ^ If might be at 
the reach of TESLA in some models with heavy neutrinos whose mass is of the order of 200 — 1000 GeV. This decay 
has also been studied within the general two Higgs Doublet Model H]. It was found that the channel Z ^ is the 
only one that may be at the reach of TESLA. Studies within the Zee Model |Q and theories with a heavy Z' boson 
with family non-universal couplings Q gave results that are far from the experimental limits. Further works have 
been realized within other models, such as supersymmetry, leptoquark theories, left-right symmetric models, etc. [|j. 
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All of the aforementioned studies have focused on specific models, which share the common feature of being of a 
weakly coupled nature, i.e. when the masses of the heavy particles become large they decouple from low energy physics. 
Therefore, it is convenient to take a more general approach that allows us to make a model-independent analysis. We 
will consider thus the effective Lagrangian approach (ELA), which is suitable for this purpose. In particular, the ELA 
has been extensively used to study some quantities that are forbidden or highly suppressed within the SM. In this 
approach there are two well-motivated schemes to parametrize virtual effects of particles lying beyond the Fermi scale 
via effective operators involving only the SM fields, namely the linear and nonlinear realizations of the electroweak 
group. 

In the linear realization or decoupling scenario it is assumed that the spontaneous symmetry breaking (SSB) of 
the electroweak group takes place in the usual way, thereby implying the existence of at least one physical Higgs boson. 
In addition, the light particles (the SM ones) fill out multiplets of SUl{2) x C/y(l). Although one can only expect 
marginal contributions from the heavy fields to low energy physics, there are indeed some processes in which the 
new physics effects may compete with the SM ones, such as those involving flavor-changing neutral current (FCNC) 
and LEV transitions |^. The latter are forbidden in the SM at any order of perturbation theory. The decoupling 
scenario is suitable to parametrize any virtual effect arising from a fundamental gauge theory that is assumed to be 
renormalizable and of a weakly coupled nature. This hypothesis is fundamental to establish a hierarchy among those 
operators of a particular dimension: gauge invariance allows us to infer the order at which the effective operators 
may be generated in perturbation theory |^. In particular, a loop-generated operator is suppressed by a factor of 
(47r)“^ with respect to a tree-level induced one. Throughout this work we will make systematic use of this fact when 
studying LEV processes mediated by the Z boson. 

As to the nonlinear realization or nondecoupling scenario fl^ , in this case it is assumed that the low-energy processes 
are affected by unknown residual strong-dynamics effects. In this effective (chiral) theory, the SSB of the electroweak 
group is accomplished by introducing a unitary matrix field U that replaces the SM doublet. It is also assumed 
that the physical Higgs boson either is very heavy or does not exist at all. The scalar sector is comprised only by 
Goldstone bosons that define the U field, which in turn transforms nonlinearly under the SUl{2) x C/f( 1) group. In 
the unitary gauge, where the Goldstone bosons are absent, we have that U = 1 and the chiral Lagrangian reproduces 
the SM without the Higgs field. Due to the fact that a strongly interacting regimen implies that loop effects can be 
as important as the tree-level ones, one can not establish a priori what operators are the most relevant. We will bear 
this fact in mind when we discuss the general structure of the Zliij couplings within the nondecoupling scenario. 

Our main aim is thus to present a model independent study of the LEV decay Z within the ELA. We 

will make general predictions for the respective rates in both the decoupling and nondecoupling scenarios. Further, 
the impact on this decay of the experimental constraints on li £jik£k and li will be analyzed, and the 

expectations at the future TESLA collider [|^, running at the Z peak (Giga-Z), will be discussed. 

The rest of our presentation is organized as follows. In Sec. || we consider the decoupling scenario and discuss 
the most general structure of the Ziitj vertex. It is argued that the contribution from the monopole structure 
7 ^ dominates over that from the dipole structure the latter can only arise at the one-loop level in any 

renormalizable theory. Therefore, the most stringent bounds (m Z ^ can be obtained from the three-body 

decay ti ijlk^k-, which receives contributions from the Zliij coupling via a virtual Z. We will also consider the 
contribution from the Zliij and WliV^. couplings to the one-loop decay li lj"f. In Sec. il, a similar analysis is 
performed within the nondecoupling scenario. We would like to stress that, in contrast to what is observed in the 
decoupling case, in the nondecoupling scenario the contributions from the monopole and dipole structures may be 
equally important due to the presence of strong-dynamics effects arising from the underlying theory, thereby allowing 
two possible scenarios. In the first case it is assumed that the monopole structure gives the dominant contribution, 
which means that the most stringent bounds on Z ^ obtained from the three-body decay li Ijlklk- 

In the second scenario it is assumed that the dipole contribution is the dominant one, which implies that, due to the 
SUl{2,) X C/y( 1) symmetry, it is possible to obtain bounds on Z ^ '^T^f by using the tree-level decays li Ij'y. It 
turns out that the bounds obtained this way are the most stringent. Finally, we present our conclusions in Sec. 


II. LFV IN THE DECOUPLING SCENARIO 

In this section we assume that the underlying theory is of a decoupled nature. The effective operators inducing 
LFV couplings were presented in a previous work 0- These operators can be classified according to whether they 
induce the ^lilj coupling or do not. 
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A. Effective operators that only induce the Ziilj vertex 


We can classify these operators in two classes. In the first place we have those operators that can be generated at 
tree-level in a fundamental theory. They are given by 



(3a) 


(3b) 


(3c) 


where Li and im stand for the left-handed doublet and the right-handed singlet of SUl{2) x respectively, 

t“ are the Pauli matrices and roman letter indices are used to denote lepton flavors. The hrst two operators induce 
the Ziiij and H£i£j couplings, whereas the third one also induces the Wiii'ij vertex. Both the Ziiij and 
couplings contribute to the one-loop induced decay £i £j'j. 

There is also another set of operators that can be generated at the one-loop level or at a higher order: 

= {UD^£rj) (4a) 


= {D^LdR,) 

Both of these sets of operators contribute to the three-body decay £i £j£k£k via a virtual Z. 


(4b) 


B. Effective operators that induce both the Zti£j and 'y£i£j vertices 


Owing to gauge invariance, operators of this kind can only arise at the one-loop level in any fundamental theory. 
According to the Lorentz structure of these operators, we can classify them in two categories: 



(5a) 

= ig' , 

(5b) 

0% = *5' {imB^^li.D,£R,) , 

(5c) 

= 9 {Ua^,W^''£R,) </), 

(6a) 

0%^ = 9' 

(6b) 


where It is understood that the hermitian conjugate of each operator is to be added in the 

respective Lagrangian, i.e. £eff. = (a*V^^)^b + H.c. We have assumed that all the effective matrices a can not be 
simultaneously diagonalized by the unitary matrices V[ and that define the mass eigenstates. Note that these 
groups of operators give rise to both Z£i£j and ^£i£j couplings as a direct consequence of the SUl{2) x Uy{£) gauge 
invariance of the effective theory. Therefore, the experimental constraints on £i —>■ £j"f can be easily translated into 
bounds on Z ^ £f£^. However, we will see below that these operators play a marginal role in this decay, though the 
situation may be different in the nondecoupling scenario. 
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C. The most general Ztilj vertex and the decay Z —> 


The effective operators shown in Eqs. (^)-(^ induce the most general Ziiij vertex. In the case of on-shell leptons, 
it is possible to make use of the Dirac equation along with the Gordon identity to transform the Lorentz structure 
induced by the operators of Eq. into a dipole structure. It turns out that the contribution from these operators 
has terms that are proportional to rriilmz or rrijlmz, with rrii^j the lepton masses. It means that these operators 
give a very suppressed contribution, as compared to that from the operators of Eq. (|^). Therefore, from now on we 
will not consider the operators of Eq. §)■ We thus can write the most general structure of the Ziiij vertex in the 
following way 




7 ^ (di ft + ft‘ift) + -J- (dift + dtft) ft + 

V / mz 7 z mz 


mz 


KPj), 


( 7 ) 


where is the Z four-momentum. We have defined the following matrices in the flavor space 


FiL = -{iy Vt (aW + a%\) V }}, (8a) 

FiR = -{yyV^aD^V^^, (8b) 

ftfi =-^^ (^) r/(c^-atw'^ + (10) 


It is evident that the terms proportional to in Eq. & do not contribute when the Z boson is on-shell. We 
thus can conclude that the contributions to Z ^ only arise from the operators given in Eqs. i) and (i)i 

i.e. only through the monopole and dipole structures. Since the monopole structure can be generated at tree-level 
by the underlying theory, its contribution will dominate that from the dipole structure because the latter can only 
arise at the one-loop level and has a suppression factor of (47r)“^. It is thus a good approximation to consider 
only the contributions arising from the operators of Eq. (^. In contrast, the ^iiij coupling is only induced by 
the operators of Eqs. (i and since the monopole contribution is forbidden because of electromagnetic gauge 
invariance, i.e. the 'yiiij coupling can only arise at the one-loop level in any renormalizable theory. In order to obtain 
bounds on Z ^ if , we will use the experimental bounds on the three-body decays ii —>■ ijikik, which may receive 
contributions from the Ziiij coupling through a virtual Z, mainly via the monopole structure. We will also calculate 
the contributions from the Zi^ij and WiiVi. couplings to the one-loop induced decay ii ij^ in order to analyze 
if this mode could be useful to obtain bounds on Z —>■ ifif- All these results can be translated readily into the 
nondecoupling scenario, where the dipole contribution to Z —> ifif will be studied also. It turns out that, in that 
scenario, the dipole contribution may be as important as that from the monopole structure. 

Taking into account just the contribution from the tree-level-generated operators, the branching fraction for the 
decay Z ^ if if can be written as 


BR{Z^ifif) 


a ('mz\ 




where we have neglected the lepton masses. We have also introduced the definition S 2 W = 2cwsw- 


( 11 ) 
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D. Bounds from the three-body decay ii —> 

The contribution from the coupling to the decay ^i ijtkik (viz Fig. |^) can be written as 

BR {I, ^ £,44) = ’ (12) 


with a = 1 — 4s^ + 8s^ and F^. being the full ii width. Again we have neglected the hnal lepton masses, i.e. 
m,j = TTik = 0. From Eqs. (0 and ( 0 ) we can obtain the following expression 



ii —!■ ijikik) 

(13) 

where BRexp.(-^z ^ ijikik) stands for the experimental constraints 



BR (/r“ ^ e“e“e“'") < 10“^^, 


(14a) 

BR (t“ ^ ijikik) < Kjk 10“®, 


(14b) 

and Kjk is a factor of order 0(1) corresponding to each different channel Q. 

These equations allows us 

to obtain the 


following bounds 


BR(Z- 

> /r=Fe±) < 1.04 X 10"^^ 

(15a) 

BR(Z- 

< 1.7 X 10"^ 

(15b) 

BR(Z - 

-4 < 1.0 X 10"4 

(15c) 


These results are in agreement with those obtained from unitarity-inspired arguments in Ref. [p^. 


E. Bounds from the two-body decay ii 

We now turn to study the contributions from the Ziiij and couplings to the one-loop decay ii ij^ (viz 

Fig. H). While the three-body decay it ijikik gets naturally suppressed by the three-body phase space and the 
exchange of a virtual Z boson, the one-loop decay ii ij^ gets a suppression factor of (47r)“^ plus an extra power 
of a. Since the current experimental constraints on both decay modes are of the same order of magnitude, the only 
way in which the radiative decay can compete with the three-body decay is if the former arises from a nondecoupling 
effect. However, we will see below that the ii ij^ amplitude is dominated by the virtual Z and vanishes when 
rriilmz 0 . 

The respective Feynman diagrams for the decay ii ijj are shown in Fig. ||. We have used the unitary gauge in our 
calculation. The expression for the Ziiij coupling was given in Eq. ( 0 )> though we will only consider the monopole 
contribution here. As for the Wiiiyi^ coupling, that can be induced by the operators of Eq. (|^, it is expressed as 

^ u{p,)PL-fi.vip,)w^, (16) 


with 


Cl = 


& 




After some calculation, the decay amplitude can be expressed as 

At'" {i^ 47 ) = u(pj) [fv - /a7®) v{pi), 

where q is the photon four-momentum and the coefficients fv,A are given as follows 


(17) 


(18) 
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fv,A = 


a 

^TTslw 




(19) 


where the superscript Vij denotes the contribution from the virtual boson {Z or W). As to the coefficients 
they are given, in terms of scalar integrals, by 


Arnfm?^ 


+ Bq — Bq — Bq) 

— 2m1 m| (±3 (i?g — iJg) + (Sg — Bg)) 


-ml {±l-Agijim%Cl) 


^IL, lfl> 


( 20 ) 


with = — 1 + 2 and 5 '^ = 2 s^. By simplicity we neglected the final lepton mass. The sign 


L (i?) term. As far as the A^^^ coefficient is concerned, we have 


(—) holds for the 




2e 


„2 

L 


[2 - 3 TO? - 2 {mly - to?) (Sg - + to? C?)] 


( 21 ) 


whereas = 0 . 


The scalar integrals Bq and Cg are given, in the notation of Ref. 


-Bo(0,to|,to|), B^ = Ro(0,TO?,to|), B^ = Rg(TO?,0,TO|), 


Bl = 


Bo{m: 


?,TO?,TO; 


.), Bl = Rg(0 


as follows: Bq = 


I ill\Y 


TOT 


^), = 


Bo{ml,0,mly), Cl = C'o(to?, 0 , 0 , to?, to|, to?) and = Co(to?, 0 , 0 , to?,/, 0 , to?^,). 

It is interesting to note that, although an effective vertex was inserted into a one-loop diagram, from the above 
expressions it is evident that the calculation renders a finite result. It can be explained from the fact that the 
coupling has a renormalizable structure. Our result is very general in the sense that it can also be applied to theories 
with an extra Z' boson with LFV couplings of the monopole-structure form. 

The branching ratio for the radiative decay li £j"/ is given by 




rrii 


Stt Ti, 


{\fv\^+\fAn = 


9 

m? a 


64 7r3 




The scalar functions involved in the coefficients A))*^ can be numerically evaluated 

We will end up with an expression of the form BR {C —> £j"f) = -\- /3i|. 

numerical coefficients. 

From the experimental side, we have 


( 22 ) 


or expanded in powers of to^. 
^ P) where the (3k are some 


BR (p - 

^ ey) < 1.2 X 10-^\ 

(23a) 

BR(t - 

-4 ey) < 2.9 x 10"®, 

(23b) 

BR(r- 

/ry) < 1.1 X 10"®. 

(23c) 


Therefore, from Eqs. (^^ and (j^) we can obtain an upper bound on the coefficients and e)(, which in turn 

can be used to put constraints on the decay Z ^ . In Fig. ^ we show the allowed region for the coefficients 

Bil ifli obtained from the decays /r ^ ey and r ^ py. As we are interested in obtaining upper bounds on 

we set e'l = 0. In the plot of Fig. the allowed regions, which interestingly are almost circular in shape, lie inside 
the curves. From these results and Eq. ( 0 ) we can obtain the following bounds 

BR(Z ^ p=^e=^) < 6.12 X I0"“, (24a) 


BR(Z^ < 2.8 X I0"k (24b) 

where £ = e ot fi. Although these bounds are weaker than the ones obtained from the three-body decay £i —>■ £j£k£k, 
they show that the one-loop decay £i —> £j^ may also be useful to obtain bounds on Z ^ £f£^. 
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III. LVF IN THE NONDECOUPLING SCENARIO 


In the scenario where the underlying new physics effects arise from a strongly interacting sector, the relevant LFV 
operators are similar to those given in the decoupling scenario, but now with the Higgs doublet replaced by the 
following unitary matrix 


/2trV“\ 

^ = “p( „ )' 

(25) 

where stands for the Goldstone bosons. In this realization of the SUr{2) x C/y(l) group, the covariant derivative 

is defined as D^[/ = i9^{7 -1- ig'Wf^U — with the Abelian field given by /2)B^. From the discussion 

presented before, it is clear that the relevant operators are the analogous of those shown in Eqs. (||) and (||), although 

(3)ii 

the operator has no nonlinear counterpart. These operators can be written as 

Cur = ^Aj)^Tr [t^U^B^U] + H.c., 

(26a) 

CuL = tX’^RTr [t^U^B^U] + H.c., 

(26b) 

Ciwu = + H.C., 

(27a) 

CiBU = g'^ {Ua^^URj) + H.C., 

(27b) 


where Ri = (0,£/{i). In this scenario there is an upper bound on the new physics scale, i.e. A ~ 4?™, which will 
adopted below. Notice that the first group of operators have dimension 4 in mass units, whereas the last ones have 
dimension 5. 

The formulas given in Eqs. ( 0 )-© also hold, but the Fi matrices are to be replaced by 

Ail = -2V[XulV^\ (28a) 

A^r = -2V^XurV^^, (28b) 

A = 2 {cw^iWU - S^XiBu) Vr ■ (28c) 


Since we are assuming a strong-interaction as responsible for the LFV effects, two scenarios are of interest. In 
the first case we will assume that the monopole contribution dominates over the dipole contribution, whereas in the 
second case we will take the dipole moment contribution as being the dominant one. 


A. Monopole dominance 


In this scenario it is 
contribution. From the 
body decay ^i 
bounds are the same as 


assumed that the structure induced by the operators (26a) and (26b) gives the 
above discussion, it is clear that the most stringent bounds can be obtained from 
. It is also clear that the relation given in Eq. (^^ still holds. Consequently, the 
those given in Eqs. Finally, the bounds arising from the decay £i —> also 


dominant 
the three- 
respective 
hold. 
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B. Dipole moment dominance 


We now neglect the monopole term and focus on the contribution arising from the operators (27a) and (27b). Due 
to SUl{‘2) X C7y(l) gauge invariance, these operators induce both the Z^i^j and vertices, which means that the 
decay £^7 can give more stringent bounds than the ones arising from the three-body decay ^ This 

is a consequence that the electromagnetic decay has a phase space factor less restricted than the three-body decay 
and does not involve the factor (rrii/mz)'^ coming from the inclusion of the virtual Z but only the kinematic one 
(jriilmzf- 

The branching fractions for the decays Z ^ and li —> can now be written as 


BR(Z 




) = 




^ ] |ylb|2 


(29) 


BR(£, ^ 



where 

^ ^ 27 rc^ + ^mu) Vr ■ 


(30) 


(31) 


After introducing the experimental constraints on the electromagnetic decays, we have 




Tz 


rui 


BRexp. (f'i —^ £j'l) 


( CW A 
V SW J 

( ®w.') 

\ cw j 


for 


‘-iwu 


for 

iui 


(32) 


By using the respective experimental constraints jj] we get 


BR{Z - 

< (8.64 X 10-2^ 7.81 x , 

(33a) 

BR{Z - 

T^e^^) < (3.09 X 10-^^ 2.79 x , 

(33b) 

BR{Z - 

r=F^±) < (1.25 X 10"^^ 1.13 X 10"^^) . 

(33c) 


where the first (second) figure in the parenthesis correspond to the operator Ciwu It should be noticed that 

the same bounds apply in the decoupling case, in the unlikely scenario where the dipole contribution dominates over 
that from the monopole. The above bounds have severe consequences. They imply that the existing experimental 
constraints on the decays ii £j^, together with SUl{2) x C/y(l) gauge invariance, are enough to rule out any 
possible detection of a LFV transition of the Z boson if it arises via a Zii£j coupling of the form of a dipole moment. 


IV. FINAL DISCUSSION 

Until now, the LFV decay Z ^ have been studied within a large variety of models iHl- These studies 
show that, at least for some values of the model parameters, the respective decay rates might be at the reach of the 
planned TESLA collider. However, all of these analyses rely on several assumptions about the parameters of the 
model under study. We have shown in this work that an ELA analysis is well suited for studying this LFV decay. 
We have considered both the linear and nonlinear realizations of the ELA. This approach has allowed us to make 
some general predictions about the Z -i- £f£^ rates starting from the current experimental bounds on the low energy 
processes U ^ £jlk£k and U ^ £jj. We summarize our results in Table 1. We also analyzed the impact of the LFV 
couplings Z£i£j and on the muon anomalous magnetic moment [g — 2 )^, but our calculation showed that 

the bounds obtained this way are rather weak. So, we refrain from showing the respective results here. In the same 
context, there are other process that could be useful to obtain constraints on the LFV Z boson couplings, such as 
/i — e conversion and muonium-antimuonium conversion. We preferred the decays £i £jlk£k and £i —> £j^ since they 
do not imply any extra assumption. 







In this work we have examined some potential scenarios for the contributions arising from the two Lorentz structures 
associated with the on-shell Ziiij vertex, namely the monopole and dipole terms. It was shown that, in the decoupling 
scenario, the decay Z ^ if arises mainly from the monopole term. In this case the strongest constraints on these 
processes are obtained from the current bounds on the decays £i —>■ £jlk£k, though the constraints on £i £jj are 
also useful for the same purpose. On the other hand, in the nondecoupling scenario, where the LFV effects have 
a strongly interacting origin, gauge invariance as the main ingredient of the effective theory induce simultaneously 
both Z£i£j and 'j£i£j vertices. In this scenario it might be that both the monopole and dipole contributions have the 
same strength. It happens that if the main contribution came from the dipole term, the current limits on the decay 
£i —>■ £jl would place severe constraints on Z ^ (^^], which clearly are far from the reach of the 

planned TESLA collider. These results suggest indeed that the dipole contribution is unlikely to be observed. 

In summary, if the new physics LFV effects are of decoupled nature, the most stringent bound on Z —)■ is of 

the order of 10“^^, which suggests that this mode would be out of the reach of TESLA. Since the current experimental 
limits on T ^ £jlk£k are less stringent than that on /i“ —> e“e“e“'", the resulting bounds onZ —> t^£'^ are also weaker 
than that on Z As a consequence, the decay Z t^£^ may still be at the reach of the TESLA collider. In 

this respect, it has been conjectured that LFV effects might be more evident in transitions involving the r lepton. 
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FIG. 1. Feynman diagram for the three-body decay it Ijlkik in the effective Lagrangian approach. The dot denotes an 
effective LFV coupling. 



FIG. 2. Feynman diagrams for the radiative decay ^ in the effective Lagrangian approach. The dot denotes an 
effective LFV coupling. There is another set of diagrams where the flavor-changing effective vertex is inserted in the opposite 
end of the Z boson or the neutrino. 



FIG. 3. Bounds on the coefficients from /i ^ ey (solid line) and t ^ fij (dashed line). The allowed region lies inside 

the curves. = 10“® and = 10“^. 
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£i ^ £j£k£k 

£i —>■ £j^ (one-loop level) 

£i ■ 

(tree-level) 

BR 

{Z^ 

7 ?^ 

i < 1.04 X 10“^^ 

< 6.12 X 10"“ 

< ( 

'10“^^ - xlO"^^) 

BR 

{Z^ 

rTe±) 

1 <1.7xl0~® 

<- 10~® 

< ( 

A0“^“ - xlO-^'^) 

BRi 

[z^ 


) < 1.0 X lO”'" 

<- 10"^ 

< ( 

'10“^^ - xl0“^^) 


TABLE I. Constraints on the LEV decays Z ^ as obtained from the experimental bounds on li ij£k£k and £i £j'y. 
The third column correspond to the monopole term of the Z£i£j coupling, which can induce the decay £i —>■ £j^ at the one-loop 
level. The last column is obtained in the scenario where the new physics LEV effects only contribute to the dipole term of 
Z£i£j. The operators that induce this term can also give rise to £i £jy at tree-level. 
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